Quasiclassical approach to high-energy QED processes in strong laser and atomic fields 



O 
(N 



A. Di Piazza^'B and A. I. Milsteini' 

^ Max-Planck-Institut fur Kernphysik, Saupfercheckuueg 1, 69117 Heidelberg, Germany 
^Budker Institute of Nuclear Physics of SB RAS, 630090 Novosibirsk, Russia 

(Dated: April 12, 2012) 

An approach, based on the use of the quasiclassical Green's function, is developed for investi- 
gating high-energy quantum electrodynamical processes in combined strong laser and atomic fields. 
Employing an operator technique, we derive the Green's function of the Dirac equation in an arbi- 
trary plane wave and a localized potential. Then, we calculate the total cross section of high-energy 
electron-positron photoproduction in an atomic field of arbitrary charge number (Bethe-Heitler pro- 
cess) in the presence of a strong laser field. It is shown that the laser field substantially modifies 
the cross section at already available incoming photon energies and laser parameters. This makes 
it feasible to observe the analogous effect in a laser field of the Landau-Pomeranchuk-Migdal effect 
for the Bethe-Heitler process. 
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The high intensity of available laser systems requires 
the investigation of the influence of laser light on funda- 
mental QED processes in an atomic field like electron- 
positron {e~^e~) photoproduction, the so-called Bethe- 
Heitler (BH) process, and bremsstrahlung This gives 
perspectives of testing QED in the presence of such in- 
tense fields, that they have to be taken into account be- 
yond of perturbation theory. The cross section of QED 
processes in a pure atomic field depends on the energy of 
the incoming particle and on the parameter Za, where 
Z is the atomic charge number and a ~ 1/137 is the 
fine-structure constant (units with h = c ^ 1 are em- 
ployed throughout). In general, the influence of the 
laser fleld is mainly characterized by the two parame- 
ters — \e\E/mujQ and x — {^/fn){E / Ec). Here, e 
and m are the electron charge and mass, respectively, 
E and Wq are the laser's electric fleld amplitude and 
photon energy, respectively, e is the energy of the in- 
coming particle (a photon in the case of the BH pro- 
cess and an electron in the case of bremsstrahlung) and 
Ec ~ vn? /\e\ — 1.3 x 10^^ V/cm is the critical electric 
fleld of QED. 

The influence of a laser field on QED processes occur- 
ring in an atomic field has already been widely studied in 
the literature, by including, however, only the laser field 
exactly (as a plane- wave), whereas the atomic field has 
been taken into account in the leading approximation in 
the parameter Za (Born approximation). For instance, 
the BH process was considered iiii2|, |3|] , bremsstrahlung 
was investigated in Q (see also the book in ffSj and 
references therein) and Delbriick scattering in f^. Since 
all above-cited results have been obtained in the Born ap- 
proximation, only the electron states exact in the plane- 
wave field (Volkov states) have been employed [1]. 

When the Dirac equation in an external field can be 
solved analytically (as in a constant field, in a plane wave, 
and in a pure Coulomb field), it is possible to obtain the 
amplitudes of a QED process via the operator technique. 



which does not imply the use of the explicit form of the 
electron's wave functions in the field. This method was 
developed in in the case of a constant homoge- 



neous electromagnetic field, in [12|, [13| in the case of a 



plane wave, and in [lj| in the case of a Coulomb field. 



When the Dirac equation cannot be solved analytically, 
the amplitude of QED processes at high-energies can be 
obtained via the operator technique based on the applica- 
bility of the quasiclassical approximation [isj. An alter- 
native method of calculation employs the quasiclassical 
Green's function of the Dirac equation in the external 
field (see [3| for the case of a pure Coulomb fieldj_ [l3 



for an arbitrary spherically symmetric field, and [18[ for 
a localized field which generally possesses no spherical 
symmetry) . 

In the present Letter we generalize the method of the 
quasiclassical Green's function to the case of QED pro- 
cesses occurring in a localized field and in a plane wave. 
We derive the quasiclassical Green's function of the Dirac 
equation, describing the propagation of an ultrarelativis- 
tic electron, almost counterpropagating with respect to 
the plane wave. This configuration is the most relevant 
from an experimental point of view. We emphasize that 
both the localized and plane-wave fields are included ex- 
actly in the calculations, which are performed in the 
leading-order with respect to the parameter m/e <^ 1. 
The obtained Green's function is then applied to calcu- 
late the total cross section of the BH process exactly in 
the parameters of the atomic and the laser field. We 
study the problem in the rest frame of the atom assum- 
ing that the energy uj of the incoming photon largely 
exceeds m, that ujqlu <^ m? (i.e., there is no pair produc- 
tion stemming from one laser photon and the incoming 
photon), and that both Za and x ~ {^/m){E / Ec) are of 
the order of unity. In this situation the total cross sec- 
tion of the pair-production process becomes independent 
of the laser parameter ^ and, since E/Ec <^ 1, the proba- 
bility of pair production only from laser and atomic field 
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|19| is exponentially suppressed. Also, we split the total 
probability of the process into two terms, the first being 
the total probability of e+e~ pair photoproduction in the 
laser field (at Za = 0), and the second representing the 
influence of the laser field on the BH process. Only the 
latter term is considered here, as the first process has al- 
ready been investigated in detail [2^, [2lj . We show that 
the presence of the laser field substantially modifies the 
cross section of BH process even at moderate values of 
the parameter x, corresponding to already available laser 
intensities and photon energies. The considered process 
indicates the feasibility of successfully applying the ob- 
tained Green's function for investigating a wide class of 
QED problems, including processes in the combined field 
of an intense laser and a highly-charged ion. 

By considering for definiteness a plane wave propagat- 
ing in the direction antiparallel to the z-axis, it is conve- 
nient to pass from the variables t and z to the variables 
(j) = t — z and T = {t + z)/2. In this way, 

"° = idt = -Pd, - Pt/2 , = -idz 



P 



'P4> + Pt/2 ., 



(1) 



D0-Pt/2, p^ 
Pt = -idr , p± = -idp , 

where p is the component of the vector r perpendicular 
to z. We first analyze the propagation of an ultrarel- 
ativistic charged scalar particle in a localized potential 
V{r) = V{p, z) and in a plane-wave field described by 
the vector potential A{t + z) — A{2T). We assume that 
the momentum of the particle is almost parallel to z. In 
this case, we can replace the coordinate z in the argument 
of the potential by T. The propagation of the particle is 
described by the Green's function of the Klein-Gordon 
equation, 



D("\T, p\T', p'- e)-- 



"('^-"^')d(°)(T, p\T', p'- e), 

^ 5{T-T')5{p-p' 



V — 



/ ds eyi]i[is{V^ -m^)\5{T ^T')5{p~ p') 
Jo 

^^'f{t, p\T', p'). 



ds ( 



= -2eV{p, T)-H, H = 2epT + [pi_ - A{T)f , (2) 

where A{T) = eA(2T). By neglecting the term V^{p, T) 
in comparison with eV{p^ T) and the term {pr, = 
Pt V -\-V Pt in comparison with epT in Eq. ([2]), we rep- 
resent the function F{T, p\ T', p') as 



F{T, p\ r, p') = L{s) e-^^"5{T ~ T')5{p ~ p') , 
L{s) 



^-ts[2eV{p^T)+H]^isH 



(3) 



The operator L{s) satisfies the equation 
dL{s) 



ds 



2ieL{s)V{Ci, T-2es), 



Ci 



2spA 



2s 



dxA{T - 2sex) . (4) 



In the leading quasiclassical approximation, we can ig- 
nore the non-commutativity of the operators V{Cij T ~ 
2es) at different values of s. As a result, we have 



L{s) — exp 



-2ise [ dxV{C2, T - 2esx) 
Jo 



C2 = P - 2sxp±_ + 2sx [ dx'A{T - 2sexx') . (5) 
Jo 



By employing the relation (see [22|) 

e-'"" = exp I - is dx[pi^ ~ A{T - 2esx)f | 

we obtain 



(6) 



-isH 



S{T - T')5{p p') = -—5{T -r- 2es) 



X exp 



Airs 



f = f dxA{T), .g2= f dxA^ir). 
Jo Jo 



(7) 



where r = T - {T - T')x. Then, we substitute Eqs. (0 
and in Eq. © and we find 

F{T, p\ T\ p') = -^J{T -T'- 2es) 



X exp 
X exp 



' 4s 



^{P-P')■f 



2ise [ dxViCi, r 
C3 = C2-2sx(/ 



X 1 



p-p' 

2s 



(8) 



where 1 stands for a constant function with unit value. 



Analogously to the case considered in [18| , it is important 
to account for the operator p± in the argument of the 
potential in Eq. ^ only for x close to xo ^T/{T - T'), 
when < < 1, otherwise the operator p^ can be 
omitted. By exploiting this fact, we can write 



exp 



-2ise [ dxViCs, t) 
Jo 



XoS 2 , 

exp [ —t- pj^ I exp 



1- Xq 

Ci= p-x{p- p') + 2sx{l - x)f , 
dy [A{Ty) - A{T'y)] . 



d^q 



2ise [ dxV{C.A, t) 
Jo 



(9) 



The useful identity 



e-'^P-ff(p) = / — ie^« g{p + 2^q), (10) 
J «7r 

valid for an arbitrary function g{p) and a positive con- 
stant 7, can be employed, with q being a vector in the 
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plane perpendicular to the vector z. In this way, the in- 
tegral over s in Eq. ([2]) can be performed, and the final 
result for the Green's function Dq{T, p\ T', p'; e) is 



the effective potential V = V + 6V, where 



i.(°)(T,p|T',p';e) = ^-^^exp 
+*(P - p')-f + i 7.^ ^ - ^rn'^so 



4so 



X / d'^qe"^ 



qe^ exp 



-2isoe / dxV{Cx, t) 
Jq 



C. = p-x{p-p') + 2^q + 2(3f, 
T 

7 , ^ = - xq)so , (11) 



T-T' 

So = 7, , 2^0 



2e 



T-T' 



with the functions /, and / being given in Eqs. ([7]) 
and ([9]) and with 9{x) being the step function. We note 
that Do{T, p\T', p'; e) does not contain anymore mo- 
mentum operators and we recall that the terms propor- 
tional to q and / are to be omitted if a;o < or xq > I. 
As it should be, the Green's function Do(T, p\ T', p'; e) 
is invariant under the replacement T ^ T' , p p' , and 
A. ^ —A.. Also, if ^ = 0, it is in agreement with the 
corresponding result obtained in the quasiclassical ap- 
proximation in 17, 3], whereas, if ^ = 0, it reduces to 
the exact Green's function in a laser field [9|. 

We consider an optical laser field (wq ^ 1 eV) and a 
particle with energy e such that ewg ^ ™^ (i-S-, e 
100 GeV). It follows from Eq. ([TTl) that in this case 
the typical value of Sq is Sq ~ l/m?, such that \T — 
T'\ ^ {e/m)\c and \p — p'\ ^ Ac, where Ac = l/m 
is the Compton wavelength. Besides, it is /3 ~ 1/m?, 
m ^ (iLoo\T ~ T'\ \A\ ^ xAc, and so{P - g^) ^ x^ 
where x = (£/™)(-E'/£'c) ^ 1 such that the effect of the 
laser field may be important. At \p + p'\ 3> Ac, we 
can neglect the terms proportional to q and / in C,x and 
perform the integral over q. The result is nothing more 
than the Green's function in the eikonal approximation. 

The Green's function of the Dirac equation can be rep- 
resented as 

G{x, x') = (x|(P _ TO + iQ)-^\x') ^{V + m)D{x, x'). 

(12) 

Here, V = ^^V^, where 7^ are the Dirac matrices, 
= p^- eA^{x), with eM' = {V{p,T),A{T)), and 
D{x, x') = {x\{V'^ - m + iQ)-^\x') is the Green's func- 
tion of the "squared" Dirac equation. It has been shown 
in [13, [i^ that it is convenient to calculate the ampli- 
tudes of various QED processes in terms of the function 
D{x, x'), and we derive here its explicit form for the field 
under consideration. We have 

V^=V^ + ia- VV{p, T) + ^i?7 • drAiT) , (13) 
where V = 9p + zdr, and >f = 7" -f 7^^. We introduce 



6V = - — 
2e 



a • VV{p, T) + ^kj ■ drAiT) 



(14) 



Although SV is small in comparison with V, it has a 
different matrix structure and the linear terms in 6V have 
also to be taken into account in the calculations (see, e.g., 
[isl). In order to obtain the corresponding correction, we 
can employ Eq. ([TT|) with the replacement V V and 
then perform the expansion with respect to SV. As a 
result we find 

D{T, p\ r, p'; e) = {1 - J-a . (dp + dp,) 

- J7 • [A{T) A{T')]}d^"\T, p\ T', p'; e) . (15) 

Note that this expression is exact with respect to the pa- 
rameters of both the laser field and the atomic potential. 

Now, we apply the Green's function in Eq. ((T5|) to 
calculate the total photoproduction cross section a in 
the laser and the atomic electric field. Due to the optical 
theorem and following Ref. [l3|, we can write a as 

cr=-Re^ de J dTdT'dpdp' 0{T -T') 

X Sp{ [(2e* • Ti^ - e*k) D{T, p\ T', p'- e)] 

X [(2e • T'l + ek) D{T', p'\T,p;e- lu)] } , 

r± = -idp - A{T) , Pi = -idp, - A{T') . (16) 

where ~ (0, e) and fc^ ~ (w, k) are the polarization 
four- vector and the four-momentum of the incoming pho- 
ton, respectively, and where the 0-function in Eq. (jlip 
is taken into account. Below, we consider the case of in- 
coming unpolarizcd photons. By substituting the Green's 
function (jisp in Eq. (jl6p and by performing the trace, 
we obtain 

a =— Re j de i dTdT'dpdp' 9{T - T') 

-{rP^Df].[V'^D^^']-^^^ (17) 
X [{Vi^ - v'l)Df] ■ [(Pi - Tl)D^°^\\ , 

where of^ ^ D^°\T, p\T' , p'; e)] and dL"^ = 
D^^\T' , p'\T, p; e — a;)]. Since we investigate the in- 
fluence of the laser field on the cross section of the BH 
process, we assume the subtraction from the integrand in 
Eq. ([T7)) of its value at A = 0. Then, we exploit the rela- 
tion cjoT ^ ujQUj/m'^ ^ 1, and expand the vector AlT) so 
that AIT) = ^(0) + £T, where the quantity € = drA 
at T = depends on some constant phase (note that 
£ = —2eE, where E is the electric field's amplitude of 
the laser wave). The final result should be averaged over 
this phase. Due to gauge invariance, the cross section a 
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is independent of ^(0). Also, it depends on E via its 
square value 5^, which in turn, for a circularly polarized 
laser field, is independent of the phase of the field. The 
integrations in Eq. (jl7|) are pe rformed in a similar way as 
described in detail in [171 , llSl ] , and we present here only 
the final result. In order to render the parameter x 
large as possible, we consider the case where uj/w? ^ rc, 
where Tc is the atomic screening radius. In the Thomas- 
Fermi model, Tc ~ [ma)^'^ Z^/^ and the cross section a 



has the form, 



cT^aa^^ix) 
a{Zaf 



ln(183Z-i/^) - ^ - f{Zc 



, f{Za) = Rc ip{l + iZa) + C , (18) 

is the Euler 



where = d\-aT{t) / dt, C = 0.577 

constant, *(x) = *i(x) + *2(x) and 



= 16Im / dx 



dre-'^f 



^!-^{x) ^ Im i dx I dre 



4 3 



-log(r)-C-*- 



19 

1 a 

15 

1 1 



,'^'^ 4 6 4/1 \ 

+ —X T a 1 - a 
15 
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-- + -a + 2ix'r'a' 1 - — a + — xV^a*(l - a) 



15 



1^223 

,ip = T+ -X a T 
32 



4a 



.16 



21 
*2(x) 



461 



7 175 



14464 
525 



xW(l-a) 



2^(1-2;), X 



15' 

uj\£\ _uj E 
2m? m Er 



Re / dy 



dx 



dre'"^ d^qln 



{q + Ff [a{2y - ifq^ + {l~a){q + Ff - ^] . (19) 



The vector F in (x) is equal to Axa\/y{\ — y)T^/'^s, 
with s being an arbitrary unit vector, = 1. Also, 
the integral over q in 5*2 (x) can be performed analyti- 
cally, but we do not report the result here for the sake 
of brevity. We note that all the Coulomb corrections 
in a are enclosed in the function f{Za), as in the case 
of a pure atomic field, and that the effects of the laser 
field are included in the functions <I>(x) and ^(x)- The 
asymptotic forms of these functions at x ^ 1 are 

^, , 28 / 2136 2\ w N 312752 , 

'^(^^ "0 + 1225^ )' *W = -^X^. (20) 

At X > 1 it is $(x) oc x^^^^ and ^'(x) oc x^^''^ Inx but 
these asymptotics are applicable only at very large val- 
ues of X- Figure 1 displays the dependence of the ratio 
cr(x)/c(0) on X for ^ = 83, where a{x) is given in Eq. 
(HH]). We have seen that this ratio depends very weakly 
on Z. The figure indicates that high-order effects in x 
become significant already at relatively small values of x- 
While the cross section results amplified for x < 0.5 (up 
to about 10% at X = 0.25), it is suppressed by about 40% 
already at x = 1 ■ This suppression is the analogous in the 
laser field of the Landau-Pomeranchuk-Migdal (LPM) ef- 
fect due to multiple scattering of a charged particle in 
matter (see the Review in [2^). Unless the LPM ef- 
fect for BH process in matter, which has never been ob- 
served because of the required ultra-high photon energies 
{u! > 2.5 TeV), the effect in the laser field can be observed 
in principle at w = 10 GeV for an already available inten- 
sity of 10^1 W/cm2 In conclusion, we have derived 
the Green's function of the Dirac equation in the qua- 
siclassical approximation by taking into account exactly 
the parameters of a strong plane wave and of a localized 
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FIG. 1. The dependence of the ratio cr(x)/o"(0) on the pa- 
rameter X for Z — 83, where cr{x) is given in Eq. (|18p . 



field. By employing this Green's function we have deter- 
mined the probability of the BH process exactly in these 
parameters. The laser field is found to suppress the cross 
section by about 40% already at x ?J 1 and this renders in 
principle feasible the observation of the analogous of the 
LPM effect in a laser field for the BH process at already 
available energy of the incoming photon and parameters 
of the laser field. Finally, the considered example shows 
the feasibility of employing the Green's function in Eq. 
(ITSI) to investigate exactly in the parameters of both the 
laser and the localized field, such processes, which so far 
have been treated in the Born approximation in Za. 
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